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Abstract. We present a constructive method to devise boundary conditions 
for solutions of second-order elliptic equations so that these solutions satisfy 
specific qualitative properties such as: (i) the norm of the gradient of one so- 
lution is bounded from below by a positive constant in the vicinity of a finite 
number of prescribed points; and (ii) the determinant of gradients of n solu- 
tions is bounded from below in the vicinity of a finite number of prescribed 
points. Such constructions find applications in recent hybrid medical imaging 
modalities. 

The methodology is based on starting from a controlled setting in which 
the constraints are satisfied and continuously modifying the coefficients in the 
second-order elliptic equation. The boundary condition is evolved by solving 
an ordinary differential equation (ODE) defined so that appropriate optimality 
conditions are satisfied. Unique continuations and standard regularity results 
for elliptic equations are used to show that the ODE admits a solution for 
sufficiently long times. 

1. Introduction 

Several recent hybrid medical imaging modalities may be recast as systems of 
nonlinear partial differential equations with known sources; see, e.g., [H |3l HI [6l 
[T6| |22| [21] for reference on such modalities. The solution of such systems requires 
that said sources satisfy specific properties which may often be recast as specific, 
qualitative properties of solutions of second-order partial differential equations. 
In the applications presented in, e.g., [HI [ID], solutions of second-order elliptic 
equations are required to have gradients that do not vanish, at least locally. In 
other applications described in, e.g., [21 [TJ [T21 UHl CH], the determinant of the 
gradients of n solutions in spatial dimension n is required to be bounded away 
from 0. 

Such qualitative properties are to be ensured by controlling the boundary con- 
ditions of the elliptic solutions. Using theories based on complex geometric optics 
solutions or on unique continuation principles and Runge approximations, it is 
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shown in, e.g., |9l HDl ES] that the quahtative properties are satisfied for an open 
set of boundary conditions that is not precisely characterized. 

This paper presents a methodology to construct boundary conditions such that 
the qualitative properties are satisfied locally. To simplify the presentation, we 
consider the setting of a second-order elliptic equation in divergence form with 
an arbitrary (elliptic) diffusion coefficient. Starting from a configuration where 
the diffusion coefficient is constant and where boundary conditions can easily be 
defined so that the qualitative property is satisfied, we propose to continuously 
deform the diffusion coefficient from the constant one to the final coefficient of 
interest. An ordinary differential equation (ODE) is then prescribed for the evo- 
lution of the boundary condition so that the qualitative property of interest is 
satisfied, at least locally in the vicinity of a finite number of points of interest, 
during the whole homotopy transformation. The qualitative property is recast as 
an adapted set of constraints. The ODE is tailored so that optimality conditions 
are met to satisfy the set of constraints. That the ODE solution exists for the 
whole duration of the homotopy transformation is guaranteed by using a unique 
continuation principle for solutions to elliptic equations. The whole procedure 
may be seen as an optimal boundary control method so that the elliptic solutions 
satisfy appropriate constraints inside the domain. 

The rest of the paper is structured as follows. The construction of boundary 
conditions ensuring that the gradient of the solution does not vanish in the vicinity 
of a given point is introduced in section [2l Section [3] presents the main results of 
this paper. In Section HJ we describe the optimality conditions that justify our 
choice of the evolution equation (the ODE) and give an example of a simpler, more 
naive, construction that does not achieve our objectives. Section |5] contains the 
proofs of the main results. Section generalizes the construction to other settings 
including the construction of solutions such that the gradients do not vanish at a 
finite number of points and the construction of solutions whose gradients form a 
basis in the vicinity of a finite number of points. 

2. Description of the Problem and Formulation of the Method 

Let X be a bounded domain in M" with boundary dX. For a given coefficient 
7(x) and a fixed x G X, the goal is to find a boundary condition / such that 
|Vti(a;)| > 1, where u is the solution of the equation 




u = f in dX. 



V ■ (7Vn) = in X 



In order to construct such an / we propose an evolution scheme. Namely, for a 
given 7o(a:) let 7^ := (1 — 5)70-1-57, Vs G [0, 1]. For a family of boundary conditions 
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{/s}se[o,i] let Us denote the corresponding solution of 

{P. 



V ■ i'jsVus) = in X 
Us = fs in dX. 



The proposed scheme consists in constructing {/s}se[o,i] with the property that 
|Vms(x)| is non-decreasing after choosing 70 and /o such that |Vmo(£)| > 1; for 
example 70 = 1 and /o(xi, X2, = Xi. Thus, / := fi is a solution of the 
problem since |Vmi(x)| > 1. 

To construct {fs}s£[o,i], let us assume that fs = fo + Jq 9tdt and denote u'^ = 
dus/ds and 7^ = d'js/ds = 7 — 70- Differentiating (P^) with respect to s gives the 
equation 

'V-(7.V<) + V-(7^Vm,) = inX 
u'g = Qs in dX. 

The condition that |Vms(x)| is non- decreasing becomes Vms(x) ■'Vu'glx) > 0. Such 
a characterization hints at the construction of {fs}se[o,i] by means of an initial 
value problem. 

We construct {fsiQs : s G [0, 1]} as the solution of 

= /o 

=0 

for an F satisfying two specific conditions. 

The first condition on the functional F is that it guarantees Vms(x) ■ V'Us(x) > 0. 
The second condition on F is that it admits a solution for initial value problem, 
with initial condition /o, for all s G [0, 1]. 

In this work, we provide an explicit description of a functional F (Definition 
13. 3p satisfying those two conditions (Theorems 13. 4[ |3I7|) . hence not only solving the 
original problem, but also providing an exphcit method to construct the solution /. 



3. Notation, Framework and Main Results 

3.1. Notation. The following notation will be used through the paper. Let X be a 
bounded domain, let dX denote its boundary, at a; G dX let v{x) denote the outer 
unit normal to X. Let X be the closure of X. The notation C'^'", /cgN, 0<a<l 
represents Holder continuity, i.e., k continuous derivatives with the k-ih derivative 
being Holder continuous of order a; in the case a=l the k-th derivative is Lipschitz 
continuous. Let C^'°'{VL) be the space of Holder continuous functions from into 
M and write X G C'^'" to mean that dX can be locally represented as the graph 
of a Holder continuous function. In C^'°'{Vt) the norm of a function / is written 
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as \ f\k,a,Q- We use the classical notation for the integrable spaces Lp{Q) and the 
norm \f\LP{n),P £ [1, oo]. Denote as W'''^{Q), k EN,p E [1, oo], the Sobolev space 
of functions with k weak derivatives in L'p{Q). In these spaces, we consider the 
usual norm that makes them Banach spaces. Let Wq '^{^) be the completion of 
C^{n) in W^'P{n); see [13j for additional details. 

3.2. Hypotheses. The following hypotheses will be assumed throughout this sec- 
tion. We assume that X is a bounded subset in M", n G N fixed. We fix p G (1, 

and let a = G (0, 1). We fix A; G N. 

We assume that X is a (7^^+^'" bounded domain, x G X is fixed. 

We assume that 7 G C'^"'"""'"^(X) and that there exist constants c,C such that 
0<c<7<CinX. 

Let 7o = 1 and /o(xi, X2, x„) = Xi. For s G [0, 1] define 7^ := (1 — 5)70 + 57. 

3.3. Main Results. The first theorem summarizes classical results and shows 
that the formal calculations in Section [2] are valid in this setting. 

Theorem 3.1. Let s ^ /, G ^^([0, 1]; C^+^'^'idX)) . For each s G [0, 1] there is a 
unique solution Us G C^+^'"(X) of the equation 

V ■ {is'^Us) = in X 
Us = fs in dX 

and s ^ Us e Ci([0,l];C^+2'°(X)). Let < = dus/dsj'^ = dus/ds and 7^ = 
d'js/ds = 7 — 70- Then u'^ satisfies the equation 




{P's) 



V ■ (7.VO + V ■ {^yus) =0 in X 
< = fs in dX 

and f (||Vu,(:r)|2) ^ yus{x) ■ V<(x). 

For y G M let (9y = (?/■ V) denote the y-directional derivative in M". Let s G [0, 1]. 
The following auxiliary problem is crucial in our analysis: 

V • (7sVA) = 9,5,. inX 
A = in dX. 



{As 



Here, 5^ is the distribution at x such that J^6xf{x)dx = f{x). The dependence 
of A on s is not written explicitly since it will be clear from the context. 

Theorem 3.2. The problem (As) above has a unique solution A G L^(X) fl 
C*^+3."(X\ {x}). IfU C (X\ {x}) IS compact, then s ^ X\u e C([0, 1]; Lp(X) n 
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We proceed to define an adequate functional F for tlie initial value problem of 
{/s}se[o,i]- 

Definition 3.3. Given f e C"=+2'"(9X) and s e [0, 1], let u e C"=+2'"(X) be the 
solution of 

'V-(7sVw) = inX 
u — f in dX. 

Let X be the solution of 

V • (TsVA) = Vu{x) ■ in X 

A = in dX. 

IfVu{x) = let > 0, otherwise let 

/_^AV-((7-7o)Vm) 
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We (Ze/ine F : C"=+2'"(9X) x [0, 1] ^ C*^+2'"(9X) as 

The functional F satisfies the required properties. 

Theorem 3.4. Given f e C^+2'°(aX) and s G [0, 1], let u e C*^+2'"(X) 6e t/ie 
solution of 

'V-(7,Vm)=0 inX 
u = f in 

Le^ = s) and let v be the solution of 

'V-(7,Vv) + V-(7;Vm) = inX 
V = g in dX. 

Then Vu{x) ■ Vv{x) > 0. 

The second property for F requires a strong relationship between the solution 
of the auxiliary problem (Ag) and its normal derivative at the boundary. In par- 
ticular, the following injectivity results is needed. 

Theorem 3.5. Let X e Lp{X) D C'=+3'"(X \ {£}) be the solution of (A,) and let 
7s9A/9z/|g^ G C'''^^'°'{dX) be its normal derivative at the boundary. Then 

dX 



A = 







dX 



y = 



6 



G. BAL^ AND M. COURDURIER^ 



This motivates us to regard A and its normal derivative as functions of y G M", 
a finite dimensional space. Using the continuous dependence of A on s, we can 
recast the previous injectivity Theorem as an apparently stronger result. 



Corollary 3.6. Let A be the solution of (Ag) and let 'ysdX/du^Q^ be its normal 
derivative at the boundary. There exists constants a,b, p,!] > independent of 
y G M" and independent of s & [0, 1], such that 



p\y\ < a 



Is 



dX 
du 



< 



k+2,a,dX 



A 



LP{X) 



< b 



Is 



du 



L^{dX) 



< v\y\ 



In particular, for any r] > 0, the quantities 

G L^{X) and 



A 






L^(dX) 



IS du 




IS du 


L'^idX) 



G C 



k+2,a 



(dX), 



as functions of y, are uniformly Lipschitz in {y G M"' : \y\ > f]}, independently of 
SG [0,1]. 

We start at s = with an adequate /o,7o, hence the estimates of the Corollary 
13.61 will imply the solvability of the initial value problem for all s G [0, 1]. 

Theorem 3.7. There exists a unique solution s ^ fg in C"*^ ([0, 1]; C'^"'"^'"(9X)) 
of the initial value problem 

'dzfs F^ifsi 



\s=0 



/o. 



In summary, for F as in Definition 13. 3[ the initial value problem admits a 
solution for s G [0, 1] and / = /i solves the original problem. 



4. Some Aspects about the Construction of F 

In this section we elaborate on the requirements on F : C^'^^'^'ldX) x [0, 1] — )■ 
C'^~^^'°'{dX) that lead us to Definition 13.31 We start by presenting a simple, naive, 
and flawed construction that exemplifies some of the difficulties before proceeding 
to the optimal aspects of Definition 13.31 

Let us consider scalings of the initial boundary condition, namely we let fs = 
0(s)/o, where G (^""^([0, 1]; M). Let Ug be the solution of 




V ■ (7s Vm,) = in X 
Us = 0(s)/o in dX. 
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Differentiating witli respect to s, u'^ lias to solve 

' V • (7s V<) + V ■ (7^ Vms) =0 in X 
w's = 0'(s)/o inSX. 

We want to construct such that V'Us(x) ■ V'U^(x) > 0,Vs G [0, 1]. Let Vs,Ws be 
the solutions of 

V ■ (7sVi;s) = in X 
Vs = 0'('S)/o in dX 

V ■ (7sVws) + V ■ iiiVus) = in X 
Ws = in dX. 

Then u'^ = Vg + Wg and (it can be checked that 0(s) 7^ 0) 

Vus{x) ■ Vm^(x) = Vus(a;) ■ Vvs{x) + Vms(x) ■ Vws{x) 

^ |Vus(£)P + Vusix) ■ Vwsix). 



and 



0(s) 

Hence, to have Vus{x) ■ Vu'g{x) > we essentially need to satisfy a condition of 
the form 

0'(s) r Vms(£) ■ Vws{x) 

— max < 0, 



0(S) I ' |Vm,(x)|2 

This condition implies the following estimate on (p 



< C\Us\k+2,a,X 



so that 

\<P'{s)\<C\<P{s)\'. 

In general, we cannot obtain any better estimate. Such an estimate guarantees 
the existence of (p for s in an open subset of [0, 1], but it does not guarantee global 
existence in [0, 1]. Indeed, the existence of for s G [0, 1] is equivalent to saying 
that for all s G [0, 1], the solution Ug of 

V ■ (7Vms) = in X 
Us = fo in dX 

satisfies Vus(x) 7^ 0. Yet, it is known that critical points of elliptic solutions do 
occur; see, e.g. [IH H [13 |2T] ■ 

This shows that may blow up in finite time if \fi.\ is large enough. We 
thus need to construct fg in such a way that \fi.\ remains sufficiently small. The 
construction of F provided in Definition l3.3l is obtained by requiring an optimality 
condition in that sense. 
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Theorem 4.1. For f G C'^+^'"(9X) and s G [0, 1] let u be the solution of 

V ■ (7,Vm) = in X 
u = f in dX. 

Let denote the solution of 

V ■ (7.Vt;3) + V ■ (7^Vm) = in X 
= g in dX. 

The construction of F{f,s) in Definition \3.3\ is such that 

F{f,s) = 8.Tgmin{\g\l,^Q^^ : g G C^^+^.a^^^^ ^ y^^^^) . > q|_ 

Proof. Let 

^ = G C'^+2'"(5X) : Vus{x) ■ Vv^ix) > 0} 
g = argmin{|^|^2(5x) : ^ ^ ^} 

From Theorems 13.41 and 13 . 7^ F{f,s) G Q, hence ^ 7^ 0. Also Q is convex and closed 
in C'^"'"^''^(9X). The objective function \g\L2(dx) is strictly convex and coercive in 
C^^^'"(9X). The existence of g does not automatically follow from this, because 
C'^"'"^'"((9X) is not reflexive, but if g exists, then it is unique. 

Theorem 13.21 and Definition 13.31 imply the existence of ft = min(0,yu) < and 
A G LP{X) C (C"=+2'"(9X))* with dX/du G C"=+2'°(9X) satisfying 

V ■ (7s VA) = Vus{x)6£ in X 
A = in dX, 

F{f, s) = fi-fsdX/du\^^ 
and 

ftVusix) ■ Vv^^^''\x) = 0. 

These are the Karush-Kuhn- Tucker (KKT) conditions |20| for the problem defining 
g. The existence of the KKT multipliers A, fi with the above conditions imply that 
g = F{f,s), and in particular imply the existence of g. The fact that the KKT 
conditions in a convex problem imply optimality is easy to check in general. We 
briefly present the calculations in this particular case for concreteness. 

If F[f, s) = 0, F{f, s) is clearly the element in Q of minimal norm. Otherwise, 
/i < and F{f, s) is such that 

Vus{x) ■Vv^^^''\x) = 0. 

We recall that for any g E Q 

Vus{x) ■ Vv%x) > 0. 
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For any g & Q, multiplication of the equation of by A and integration by parts, 
gives 



dX 

idx du 



Vu^ix) ■ Wix) = / AV ■ (7:Vn,) - / 
Subtracting this last expression for g & Q and F{f, s), we obtain 
- / 'y^ira + / 7sTr^(/, ^) = VM,(a;) ■ Vv\x) - Vusix) ■ Vv^^f''\x) > 

JdX OP Jqx ou 

f dX^,^, f dX 
JdX du Jqx du 

Since fi < and F{f, s) = fi'jsdX/di' ^ the previous inequality implies 

dX 2 



76 



L^{dX) 



< 



Is 



dX , 

dx du 
dX 



L2{dX) 



L^{dX) 



dX 
du 



L2{dX) 



< \g\L^(dX) 



^ \Fif,s)\L2(ax) < \g\L'^(dx) 
proving that F{f,s) is the element in Q of minimal L'^{dX) norm. 



□ 



In summary, among all the possible choices of F satisfying the non-decreasing 
norm of the gradient at x, our definition of F{f, s) is the one of minimal L'^{dX) 
norm at each s G [0, 1]. 

5. Proofs and Intermediate Results 

5.1. Proof of Theorem 13.11 In this subsection, let A; G N fixed, < a < 1 fixed. 

Theorem 5.1. Let X be a C'^+a." bounded domain m W. Let f G C'^^'^idX) 
and h G C^'"(X). Let 7 G C^^^'"(X) be such that 3c, C constants for which 

< c < 7(x) <C <oo Vx G X. 

Then there is a unique solution u G C'^+^'"(X) of the equation 

V ■ (7VM) = h in X 
u = f in dX, 

andu satisfies the following estimate where the constant k, depends only on n, a, c, C 
and X , 

+ \h\ 
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For the previous Theorem, existence is estabhshed in [T3l Thm. 6.14, Thm. 6.19, 
Lem. 6.38] and the estimate is a consequence of \13\ Thm. 6.6, Lem. 6.38, Thm. 
3.7]. The estimate and the hnearity of the problem imply a smooth dependence 
of u with respect to the boundary condition and the equation coefficient. This is 
stated explicitly as follows. 

Corollary 5.2. Let X be a (7'^+^'" bounded domain in M". For an interval / C M 
let s ^ fs e _C'i(/;C"=+2'°(9X)) and s ^ hs e Ci(/; C"='"(X)). Let s ^ 7, G 
C^{I] C''^^'"'{X)) and such that 3c, C constants for which 

< c < -fs{x) < C < 00 Vx G X, Vs G /. 

// we let Us, s G /, be the solutions of 

V ■ (7sVms) = hs in X 
Us{x) = fs{x) in dX, 

then s ^ Us e C^I^C'+^X)). Letting 7^ := d-fs/ds , f'^ := dfs/ds and h'^ : = 
dhs/ds, we also get that u'^ := dus/ds satisfies the equation 

V ■ (7sV<(x)) + V ■ {iyus{x)) = K in X 
Kix) = f'six) in dX. 

In addition, for a given x E X , we have ^^|VMs(x)pj = 2Vus{x) ■ Vu'^^x). 
Proof of Theorem 13.11 It is a direct consequence of Theorems 15 . 1 1 and Corollary 

□ 

Remark 5.3. Theorem 15.11 and Corollary 15.21 remain true if we replace V ■ (7V) 
by any uniformly elliptic operator L = aijOx^d^^ + hid^^ with ajj, hi G C^'°'{X). 

5.2. Proof of Theorem 13.21 In this subsection let A; G N fixed. Let p G (1, ^) 
fixed. Let a = (n - |) G (0, 1). 

For ?/ G M let 9y = (?/ ■ V), s G [0, 1], we study the auxiliary problem. 

\ .{j,V\) = dy5i inX 



{As 



A = in dX. 



Intuitively, the solution A is a directional derivative of a Green's function, and so it 
should behave as a Green's function with one degree less of regularity. Among the 
statements in Theorem 13. 2[ the uniqueness of A is the simplest and follows from 
standard arguments. The continuous dependence of A on s is the most technical 
aspect and it will require the explicit construction of the singular part of A. This 
construction will also prove the existence and regularity stated in Theorem 13.21 
The construction of the singular part of A is presented in a couple of technical 
lemmas below. We start by introducing the necessary notation. 
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Definition 5.4. Let E G 'N. We say that {cj}j(zE is a family of homogeneous poly- 
nomials centered at x if each Cj is a polynomial formed exclusively by monomials 
centered at x of total degree j , namely 

m=j 

where (3 G W, = Y.l=iPi> - = ^1=i{xi - ^if' and each c^j G M. We 
say that {c^,j}|/3|=j C M are the (finitely many) coefficients of cj. 

Definition 5.5. Let E = N or E = {0, 1, 2, ra}. Let {cj}j^E be a family of 
homogeneous polynomials centered at x with cq ^ 0. We define the family of 
functions {vm}m(^E associated to {cj}j^E as follows. 

Let B be an open ball centered in x and containing X. Let g be the solution of 



Ag = Sx in B 
g = in dB. 



Then define vq := -^dyg. 



Let w be the solution of 



Aw = Vq in B 
w = in dB. 



Then define Vi := ^(Vci ■ Vw — CiVq). 



For 2 < m,m G E , define recursively Vm o-s the solution of 

H=0 



CoAvm = J2T=o^['^ ■ (viVcm-i) - A{cm-iVi)] in B 
Vm = in dB. 



Lemma 5.6. A family {vm}meE from Definition \5.5\ satisfies 

(a) voe LP{B)r]C^{B\{x}). 

(b) djVo G W^'^{B) for any dj homogeneous polynomial of degree j centered at 
X, Vj G N. 

(c) vx G n C°°(5 \ {x}). 

(d) We get in X 

V ■ {coVvo) = dy5i (5.1) 

m 

^ V ■ {cm-iVvi) = 0, Vm > 1 (5.2) 

1=0 
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Proof. Using the same notation as in Definition 15. 5[ g is the Green's function 
at X of the Direchlet problem for the Laplacian in B, hence there is an exphcit 
expression of g; see [13]. Since Cq is constant and Vq := -^dyg, properties (a) 
and (b) are automatically verified from the explicit expression for Vq (recall that 

Property (a) and the definition of w imply w G W'^'^^B) fl C°°{B \ {x}); see [T3| 
Thm 9.15]. Then (c) follows from the definition of Ui and (b). 

Property (d) is the definition of rewritten. □ 

The family {fmjmGE has the following regularity. 

Lemma 5.7. Let {cj}j(zE, {dj}j(zE be families of homogeneous polynomials cen- 
tered at X. Let {fmjmeE be the family of functions associated to {cj}j(zE in Defi- 
nition 15.51 Then 

(1) Vm e VV^'PiB) n C^{B \ {£}), Vm E E. 



(2) djV^ e Vr^+'"'P(5),Vm e E,Vj > 1. 

(3) dj/^Vm e Vr™+^-2'P(5),Vj > l,m > 1. 

Proof. The proof is by induction. As the base case, from the previous Lemma we 
already have (1) for m = and (2) for m = 0, Vj > 1. We also have (1) for m = 1. 
The following steps complete the induction argument. 



(1)V0 < m < M and (2)Vj > 1,V0 < m < M 
Using the definition of vm, M > 1 

d ^ ^ ""^ 

djAvM =— ^ [V ■ (viVcM-i) - A{cM~-iVi)] 



(3) for m = M, Vj > 1 



Co 



i=0 
M-1 



1 V— \ 

'~ ■ (djViVcM-i) - ViV ■ {djVcM-i) - djA{cM-iVi)] 

Co ~: 



1=0 
M-1 



= — 52 [V ■ {djVi\/cM-i) - ViV ■ {djVcM-i) - A{djCM-iVi) 

^0 i=0 

- CM-iViAdj + 2V ■ {cM-iViVdj)] 
and by the inductions hypotheses each summand in the right hand side is in 

^i+M-2,p(^)^ 



For M > 1, 



(1) and (3) for j 



(2) for J = 1 



We have 



AdiVM = 2V ■ {vMVdi) + diAvM 
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and by induction hypotheses the right hand side is in W^^ ^'^{B), hence (Chp. 9, 

[15]) diVM e W^'+^'P{B). 



For M > 1, 



(2) for J 



We have 



(3) for J = J and (2) for 1 < j < J = 

AdjVM = 2V ■ {vM^dj) + djAvM - VM^dj 

and by induction hypotheses the right hand side is in hence (Chp. 

9, [15]) djVM e W^+^'P{B). 



(1)V1 < m < M and (2)Vj > 1,V1 < m < M 
definition 

AI-l 

Avm = — X] [V ■ (ViVcM^i) - A{cM-iVi)]. 



1) for m = M 



• By 



Co 



i=0 



By induction hypotheses, forM > 2 the right hand side is in W'^^ ^'^(5), hence 
Vm £ ly^^'^(-B). Elhptic regularity and the induction hypotheses also imply Vm ^ 
C°°{B\{x}). □ 

In Definition 15.51 we have an explicit construction of each Vm in terms of the 
polynomials Cj. This provides an explicit dependence of each Vm in terms of the 
coefficients of the Cj's. 

Lemma 5.8. Let {Vm}meE be the family associated to {cj}j(zE- We can write each 
Vm as 

where Im is a finite index set, {pi,m\ieim ^■^ a family of real valued polynomials eval- 
uated in {I/cq} U {c/3j}|/3|=j,i<j<mj but otherwise independent of {cj}. And where 
{cz,m}/e/m ^ family of functions in X independent of {cj}, each ci^m satisfying 
(1),(2),(3) for m of Lemma\^ 

Proof. By induction. True for Vq from its definition with pi o(l/co) = I/cq and 
Ci,o = dyQ. For m > 1, the linear system defining Vm can be written as 

^ m—1 

Avm = V(cm-iVt;i) 



^0 7=0 

m—1 



AVr. 



Yl C^,m-Niix-xfVl 



Co t 



i=0 \l3\=m-i 



m—1 



E E E 

i=0 |^|=m-i leh 



Cis- 



co 



'-pi,,Viix-xfVe^). 
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Defining {ei^rn}ieim tlie solutions e of tlie family of equations 

Ae = V((a; - x)'^Vei^i) foT < i < m - = i,l e li 
the result follows. □ 
We can now explicitly describe the singular part of the solution A of (Ag). 

Theorem 5.9. Let 7 G C^{B). Let {cj}f^Q form the partial Taylor sum of 'j 
about X, namely 

K 
j=0 

with 'Jk{x) = o{\x — x\^),^K ^ C^{B). Assume Cq ^ and let {Vm}m=i 
the family constructed in Definition \5.5[ corresponding to the {cj}f^Q. Define 
wk e n C^{B \ {x}) as 

K 

= ^ Vm- 

r?i=0 

Then there exists Hk G W^-^'P{B) n \ {x}) such that 

V ■ (7VWA') = dy^x + hx in B. 

In addition, if U is a compact subset of B \ {x}, then Wk G L'^{B) fl C°°{U) 
depends continuously in the coefficients of {c^jj^Q. ^4/50, hx G fl 
C°°{U) depends continuously in 7 under {B) perturbations. 



Proof. Let 

K-i 
j=0 

Then jK-i{x) = o{\x - x\'^-'),-fK-i e C^(B). We have 

K 

m=0 
K K-m 

= X V ■ Q X Cj- + -fK-m VW„ 
m=0 j=0 
K K-m K 

= X X V ■ {CjVVm) + X] ^ ■ ilK-nyVm) 
m=0 j=0 m=0 
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K i K 

1=0 j=0 m=0 

= dyS^ + hK 

where the first term is simphfied using equations (15. ip . (15. 2p and where hx is 
defined as 

K 

m=0 

Then G W^~'^'^{B)r[C°°{B\{x}) by Lemma IFTTl The continuous dependencies 
of wk and Hk are a consequence of Lemma [5.81 and the definitions of wk-, hx- D 

Theorem 5.10. Let X be a C'=+2-" bounded domain m W. Let 7 G C^+"+2(X) 
be such that 3c, C constants for which 

< c < 7(x) <C <oo Vx G X. 
Then there is a solution X G Lp{X) f] C^+^'^^iX \ {x}) of 



(A) 



V ■ (7VA) = dy6£ in X 
A = in dX. 



Also, for any compact set U C {X \ {x}), we have that X\u G Lp{X) fl C*=+^'"(f/) 
depends continuously in 7 under C*''+"+^(X) perturbations. 

Proof. Let i? be a ball centered in x and large enough to contain X. Extend 7 as 
(jk+n+2^j^^ and let if = A; + n + 2. Let {cjjj^o form the partial Taylor series of 7 
about X and let wx, be as in Theorem 15.91 Since Hk G then (Chp. 

9, |15j) there exists a unique v G l^^'^(-B) solution of 

V ■ (7VtO = -/^A- in 5 
f (x) = in dB 

which depends continuously on Hk- By Sobolev embedding, v G C''~^'^'°'{B) (re- 
call that a = n — n/p) and it depends continuously on Hk, hence it depends 
continuously on 7 under C^~^^~^'^{B) perturbations. 

Additionally, since [wk +v]sx ^ C'^"'"^'"(9X), Theorem 15. 11 implies that there is 
a unique w G C^+^'"(X) solution of 

V ■ (7Vw) =0 in X 

w{x) = —wk{x) — v{x) in dX 

which depends continuously on [wk + v\qx ^ C'^^^'°((9X), hence it depends con- 
tinuously in 7 under C^^"'^'^{B) perturbations. 

Finally, A = ^ solution of {A) with the desired properties. □ 
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Proof of Theorem 13. 2L Theorem 13.21 is Theorem 15 .101 for k + 1 instead of k. □ 

5.3. Proof of Theorem 13.41 We use the same notation as Definition 13.31 and 
Theorem 13.41 

Proof. We separate in two cases. 
Case 1. If Vu{x) = then immediately Vu{x) ■ Vf (x) = 0. 

Case 2. The equations for u,v and A, plus integration by parts, give 



Case 2.1. If > then 51 = and AV ■ ((7 - 7o)Vm) > 0, hence 
Vm(x) ■ Vf (x) = AV ■ ((7 - 7o)Vm) > 0. 

Case 2.2. If < then g = ndX/du and we get Vu{x) ■ Vf (x) = 0. 

□ 

5.4. Proof of Theorem 13.51 We use the notation of Theorem 13.51 

Proof. It is clear that [?/ = O] ^ [A = O] ^ ['jsdX/du^Q^ = O] . In the opposite 
direction. Assume 7s(9A/9z/L„ = 0, then A satisfies the equation 



By unique continuation A = in X \ {x}. Since A G U'{x) we conclude A = 0. 
Finally, liy^Qlei^e C^{X) be such that dy^{x) ^ 0. Then A ^ G U'{X) 



5.5. Proof of Corollary 13.61 We start with a lemma about injective linear maps 
defined over a finite dimensional domain. 

Lemma 5.11. Let I dM. he a closed hounded interval. Let (V, | ■ [y) he a normed 
vector space. Let if^ : M" — i- V, s G /, he a family of injective linear functionals. 
Assume 






V ■ (7,VA) = in X \ {x} 
A = mdX 
7s 1^ = in dX. 



since AV ■ (7^ Vy^) = dyLp{x) ^ 0. 



□ 



lim Ht{y) = H,{y),yyeW,\/seI. 
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Then there exist constants < a,b < oo such that Vs G / 

a\y\<\Hsy\v <b\y\, Vy G M". 

Proof. Let {cj}"^]^ be a basis of M". Since / 9 s i— )■ Hs{ei) are continuous and / is 
compact, maxj=i ... „ sup^gj iHgCil < oo. Since Hg are linear, the existence of 6 > 
for the second inequahty follows. 

Assume $a > such that the first inequality holds. By the compactness of 
/ and the linearity of each Ht, there exists s G / and I 3 t s, together with 
yt G M", \yt\ = l,yt ^ ys, such that \Htyt\ ^ 0. But then (using \Ht{ys-yt)\v < 

b\ys-yt\) 

< \H,y,\v < \H,{y.,) - H^lv + \Ht{y., - yt)\v + W^y^W ^ 0. 

Hence Hgys = 0, contradicting the injectivity of Hs since \ya\ = 1. □ 

Proof of Corollary 13.61 From Theorem 13.51 the linear maps M" 9 y i— i- A G 
LP{X) and W 3 y ^ -fsdX/du G C"=+2'°(9X) C L'^{dX) are injective Vs G [0, 1] 
(A is the solution of (As)). From Theorem 13.21 for ?/ G M" fixed, A G U'{X) and 
'jsOX/du G C^+^'°(9X) C L'^{dX) depend continuously on s G [0, 1]. Lemma [5.111 
then implies that all the quantities 

\y\, \Mlp{x), \lsdX/du\k+2,a,dx and \'jsdX/du\L2(gx) 

are comparable uniformly Vs G [0, 1]. The last statement of Corollary l3.6l is true for 
any quotient of two Lipschitz function in a set where the denominator is bounded 
away from zero. □ 

5.6. Proof of Theorem [Ml Let us recall the definition of F : C"=+2'"((9X) x 
[0, 1] ^ C+^'^'idX). Given / G C^+2'"(aX) and s G [0, 1], let u G C"=+2'°(X) be 
the solution of 

V ■ (7,Vm) =0 in X 
u = f in dX. 

Let A be the solution of 

V ■ (7sVA) = Vm(x) ■ in X 

A = in dX. 

If Vu(a;) = let /i > 0, otherwise let 

/^AV-((7-7o)Vm) 



/i - 2 



L2(9X) 
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We defined 

if /i > 
if < 0. 

Lemma 5.12. There exists a constant k > independent of s E [0, 1] such that 

k+2,a,dX < l^\f\k+2,a,dX- 

Proof. 

When F{f, s) = there is nothing to prove. Otherwise (A ^ 0) 

d\ 

k+2,a,dX 



□ 



|-^(/) s)\k+2,a,dX 



LAV 


((7- 


- 7o)Vm) 






2 

L2(9X) 



' du 



k+2,a,dX 



< 



< 



X 



X llsQ^lL^idX) 



V-((7-7o)Vn) 



|7sa;;|A:+2,a,ax 



nsQ^iL^idx) 





A 


LP(X) 


7 ^1 

Is 


k+2,a,dX 


1 '■•^du 




1 lsd„ 


ll^idX) 



V- ((7-7o)Vm) 



Lp/(p-i)(X) 



We used Holder inequahty to go from the third to the fourth hne. Corollary 13.61 
and the boundedness of X to go from the fourth to the fifth line, and Theorem 
15.11 to go from the fifth to the last line. 

Definition 5.13. Given r]>0 ands e [0, 1] let us define the set Nr,,s C C^+2.°(5X) 
as follows, f G Nfj^s if o-nd only if the solution u of the equation 

V • (7,Vm) = in X 
u = f in dX 

satisfies |V'u(x)| > t]. 

Lemma 5.14. Fix i] > 0, for f e Nr^^s C C''+^'°(9X) and s E [0, 1] let u, A and /i 
be the ones involved in the definition of F{f,s). Then 

• Nri^s 3 / >■ ^/\ls^\L^(dx) £ L^{X) is Lipschitz continuous and bounded, 
uniformly in s G [0, 1]. 



^»7,s 3 / 7sf~/l7sf~|L2(gx) £ C is Lipschitz continuous and 

bounded, uniformly in s E [0, 1]. 
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• Nfj^s 9 / I— )■ M G C'^~''^'°(X) is linear continuous, uniformly in s & [0, 1]. 

Proof. The last property is a direct consequence of Theorem 15.11 The first two 
properties are quickly deduced from Theorem 15. 11 the definition of A and Corollary 
13:61 □ 



Theorem 5.15. Given r] > there exists k > such that^s E [0, 1], V/i, /2 G A', 



r],s 



S) - F(/2, s)\k+2,a,dX < + \ fl\k+2,a,dX + \ f2 



k+2,a,dX 



)l/l-/2| k+2,a,dX- 



Proof. Let Ui, Aj, /ij, i = 1, 2 be the values appearing in the definitions of s) 
and F(/2, s) correspondingly. 

If /xi, /i2 > then F{f\, s) = F{f2, s) = and \F{f\, s) - F{f2, s)\k+2,a,dx = 0. 

If /ii > and /i2 < then s) = and 



- F{f2,s)\k+2,a,dX = \F{f2,s)\k+2,a,dX 

\2 



X lls^ll^dX) 



V- ((7-7o)Vm2) 



\'ys^\k+2,a,dX 



llsf^ll^dX) 



<P 



A, 



X hs^lL^idX) 
A2 

X lls^ll^dX) 



-V ■ ((7 - 7o)Vm2) 



V ■ ((7 - 7o)Vm2) 



Ai 



1^, I 
X ns-Q^\L^{dX) 



V-((7-7o)Vni; 



From the second to the third line we used Corolarry 13. 6[ From the third to the 
last line we used the fact that each integral has the same sign as the corresponding 
yUj, and we are in the case of HiS with opposite signs. 
If /ij, /i2 < then 



\F{f^,s)-F{f2,s)\k+2,a,^X 



Ai 



X ns-Q^\L'2{dX) 



V- ((7-7o)Vui) 



du 



U, dXi I 
Hs^lL^idX) 



A, 



X lls^ll^dX) 



V- ((7-7o)VM2)- 



^ ^ 



hs^lL^dX) 



k+2,a,dX 



Using Lemma 15.141 we observe that in any of the three cases, we are left with 
products of bounded Lipschitz functions and one continuous linear function, all 
bounds being uniform in s G [0, 1], which readily implies the estimate above. □ 

Proof of Theorem 13.71 Let > be such that /o G A^,,,o and let p > be such 
that \F{f,s)\k+2,dx < p\f\k+2,dx for all / G C''+\dX),Ws G [0,1] (such p exists 
by Lemma [5.12p . 
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In order to prove Theorem 13.71 it is enough to show that there is {/s}se[o,i] ^ 
C([0, 1]; C"=+2,«(ax)), such that Vs G [0, 1], 

/. = /o+ f F{U,T)dT. (5.3) 
Jo 

Writing the initial value problem in this integral form, the uniqueness of the solu- 
tion will be consequence of the existence proof (Step 2 below, which uses a Banach 
fixed point argument), the continuous differentiability in s will be automatic from 
the continuity of s t— )■ and the continuity of F (Theorem 15. ISp . 

To prove that there exists {/s}se[o,i] satisfying Equation (15. 3p Vs G [0, 1], we fol- 
low the proof of Picard-Lindelof Theorem for ODEs with some small modifications. 
The proof is done in two steps. 

Lemma 5.16 (Step 1). Let Q < t < I. If {fs}sem C C^+'^'''{dX) satisfies 
Equation / fO|) Vs G [0,t), then ft := {Ymvs^t- fs) G C^+^'"(i9X) exists and ft G 
Nj^^t (starting with /o G N^^fl)- 

Proof of Step i. If Vs G [0,t) 

fs = fo+ f F{f,,T)dT, 

Jo 

then Vs G [0, t) 

\fs\k+2,a,dX < \ fo\k+2,a,dX + / |-^(/r, t) |fc+2,a, 

dxdr 

Jo 

^ \fo\k+2,a,dX + P / | /r | fc+2,a,aX'^T 

Jo 

hence \fs\k+2,a,dx < e'"'|/o|fc+2,a,ax, Vs G [0,t). In particular, for < Si < Ss < t 

\fs2 fsi \k+2,a,dX < / \F{fT,r)\k+2,a,dxdT 
J si 

< pe''*|/o|fc+2,a,ax|s2 - Si I, 

i.e., {fs}se[o,t) is a Cauchy limit as s — > . Since C^~'"^'"(9X) is complete, ft := 
(lim^^t- fs) G C''^^'°'{dX) exists. The inequality for \fs2-fsAk+2,a,dx also implies 
{fs}se[o,t] continuous on s G [0,t], hence continuously differentiable in s, and since 
/o G N^fl, Theorem [331 implies fs e iV^,s,Vs G [0,t]. □ 

Lemma 5.17 (Step 2). /// G N^^^t then there exists e > and a unique {fs}se[t,t+e) £ 
C{[t,t + e)] C"=+2'"(aX)) s^ic/i that 



f^ = f + I nfr.r)dr, VsG [t,t + e) 
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Proof of Step 2. Let us recall that / G N^,t C C^+'^'''{dX) if and only if the 
solution u of the equation 

V ■ (7tVM) = in X 

u = f in dX 

satisfies |Vn(i:)| > r]. The smooth dependence of u in terms of the boundary 
condition and the equation coefficient (Theorem IS.ip . implies the existence of 
e > and 5 > such that if /i G C^+^'"((9X) satisfies \h\k+a,ax < ^, then 
+ G iV,,„Vs G [t,t + e). 

Let us consider the following; non- empty closed set of C([t, t + e); C'=+2."(5X)) 
^={{/sKeM+.)eC'([t,t + e);C'=+2,"(ax)) :/i = /, sup \fs - f\k+2,a,dx < 5}. 

s£[t,t+e) 

If e > is small enough, we can define the following operator T : T ^ T 

T({/,}). := r({/,},e[M+e)). '"^^ e ^ + 

Let us verify that {T({/,})4 G ^ for {/,} G ^. First T{{fs})t = /, also 



< tp sup \fr\k+2,a,dX 

< ep(|/|fc+2,a,ax + 5) 

< 5 (if < e small enough) , 

and for t<(j<r<t + e 

\T{{fs})a -T{{fs})r\k+2,a,dX < \F {fr, r)\k+2,a,dxdT 

<W- r\p{\f\k+2,a,dX + S) 



and hence a ^ T({/,})^ G C"=+2-"(5x) is continuous and {T({/,})^} G J". 
In addition, if {/«}, {(7s} G J-", then for any a G [t, t + e) 

l^({/J)<T - T({5(J)^|fc+2,a,ax < ^ \F{fT,r) - F{g^,T)\k+2,a,dxdT 

< eL sup [(1 + \ fs\k+2,a,dX + |5's|fc+2,a,ax)|/s " 5's|fc+2,a 

se[t,t+e) 

< efi;(l + 2|/|fc+2,a,ax + 2(5)) sup \fs - gs\k+2,a,dx 

s£[t,t+e) 

< sup I/, - gs\k+2,a,dx (if < e small enough). 

^ s£[t,t+e) 
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To go from the second to the third hne we used Theorem 15.151 and the fact that 
fs, Qs e N^^s for all s e [t,t + e). 

Hence, we have that T : J-" — i- J-" is a contraction in a closed subset of a Banach 
space. The Banach fixed point Theorem implies the existence of a unique fixed 
point, hence the proof of Step 2 is complete. □ 

Putting together Step 1 and Step 2, since [0, 1] is connected, we conclude the 
existence of a unique {/s}se[o,i] that solves Equation (15. 3 p for all s G [0, 1] (start- 
ing with /o G A^r?,o for some > 0). This family is continuous in s, therefore 
continuously differentiable in s. Completing the proof of Theorem 13.71 □ 

6. Extensions 

The evolution scheme presented in the previous sections solves constructively 
the following problem: given a smooth enough bounded domain X and coefficient 
7, and given any point x G X, find a boundary condition / such that the solution 
u of 

fv ■ (7Vm) =0 in X 
\u = f in dX 

satisfies |V'u(x)| > 1. We now consider two possible extensions, one that imposes 
a condition over finitely many points instead of only one, and one that imposes a 
condition involving finitely many equations. 

6.1. Finitely Many Points. Given a bounded domain X, a coefficient 7 and 
finitely many different points {xijjg/ C X, the goal is to find a boundary condition 
/, such that the solution u of the equation 

V ■ (7V'u) = in X 
u = f in dX 

satisfies |Vu(a;j)| > l,Vz G /. The process is analogous to the case of one point. 
We are now considering multiple constraints to be satisfied although we have only 
one equation and one boundary condition to control. The scheme proposes to start 
with an appropriate 70, /o (e.g. 7o = 1 and /o(xi, X2, a;„) = xi) and construct 
fs such that the solution Us of 

V ■ (7.Vms) = in X 
Us = fs in dX 

and the solution u'^ of 

V ■ (7s V<) + V ■ i-f'sVus) =0 in X 
< = fs in dX 
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satisfy ^|VMs(a;j)p = 2Vus{xi) ■ Vu'^^Xi) > 0,Vi G /. Again, we construct fg as 
the solution of an initial value problem 

fs\s=0 = fo 

for an appropriate F defined below. 

We construct the functional F : C"=+2'°(5X) x [0, 1] C"=+2-"(5X) as follows. 
We assume X is a (7^+^'° bounded domain. Let 70,7 G and let 7^ = 

[(1 — s)7o + S7], s G [0, 1]. Assume < c < 70, 7 < C < 00. For s G [0, 1] and for 
/ G C*--+2.°(5X) let u be the C^'+2.°(X) solution of 

fV-(7sVM) = inX 
\u = f in dX. 

Let V be the C'=+2'"(X) solution of 

V ■ (7,Vt;) + V ■ ((7 - 7o)Vm) = in X 

V = g in 9X 

where (7 G C'^^^'"(9X) will be prescribed below. 

The difference with the previous process appears in that we need to consider 
many auxiliary problems. Let Xi,i G /, be the solutions of 

V ■ (7.VAi) = Vu{xi) ■ V(5£, in X 
Ai = in dX. 

From Theorem [321 k G Lp{X) and {-fsdXi/du) G C^+2.°(5X) depend continu- 
ously on s. Since the are different, the {Aj}jg/ are linearly independent 
and the {7s9Aj/9z/}jg/ are linearly independent (as long as Vu{xi) 7^ 0,Vi G /). 
This is proved exactly as in Theorem 13.51 

By integration by parts, we obtain for alH G / that 

Vu{xi) ■ Vv{xi) = f XiV- ((7 - 7o)Vn) - / Is^g, 
Jx JdX du 

and we let g G span({7s9Aj/9z/}jg7) C C''"^^'"(9X) be such that, e I 

f ^^ir9= f A,V-((7-7o)Vn). 
JdX Jx 

Hence, Vu{xi) ■ Vf (xj) = 0, Vi G /. Since the {'jsdXi/dh'}i^i are linearly indepen- 
dent in L'^{dX), such a g exists and is unique (as long as Vu{xi) 7^ 0, Vi G /). 

Also, by an extension of the finite dimensional argument leading to Corollary 
13.61 (adding the linear independence of the {■jsdXi/di'}i(zi), there exist constants 
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Ci, C2, independent of s and /, such that 



<Ci 

k+2,a,dX 



U 



<C2 

2,X 



f 



k+2,a,dX 



By defining F{f,s) := g, the continuity and boundedness of F : C''^^'"{dX) x 
[0, 1] — > C'^+^'°(9X) can be proven exactly as it was done in the previous case 
(observing that the evolution with F keeps |Viis(xi)| > l,Vi e /). This provides 
the following result: 

Theorem 6.1. Assume X is a C'^+^'°' bounded domain. Let 70,7 G C'^+"'+^(X) 
and let 7^ = [(1 — 5)70 + sj],s G [0, 1]. Assume < c < 70,7 < C < 00. Let 
/o,7o be chosen appropriately (e.g. 70 = 1 and fQ{xi,X2, ...,Xn) = xi) . Define 
F : C'^+^'"(9X) X [0, 1] C'^+^'"(9X) as above. Then there exists a unique 
solution {/s}se[o,i] ^'■'^([0, 1]; C'^+^'"(9X)) of the initial value problem 

gj/s — Fifs-i 

/s|s=0 = /o- 

The family {fs}s£[o,i] constructed in this way, satisfies that each Ug, solution of 
(Pg) with boundary condition fs, is such that \Vus{xi)\ > 1, Vi e /, Vs e [0, 1]. 

Hence / = fs\s=i solves the problem presented at the beginning of this Subsec- 
tion, with a condition imposed over finitely many points. 

6.2. System of Equations. Given a bounded domain X C R^, a coefficient 7 
and fixed point a; G X, the objective is to find boundary conditions {/*}i=i,2,3, 
such that the solutions k\ i = 1,2,3 of the equations 

'V-(7ViiO = inX 

u\x) = r{x) in dX 

satisfy det[(VM*(a;))j=i,2,3] > 1- We now have only one constraint to satisfy. It 
involves multiple equations and multiple boundary conditions. The scheme pro- 
poses to start with appropriate 70, {/o}i=i,2,3 (e.g. 70 = 1 and /o(xi, 3:2, X3) = xi) 
and construct /], i = 1, 2, 3 such that the solutions u\ of 

V-(7.V<) = inX 

ui{x)^fi{x) indX 

and the solution (ul)' of 

V ■ (7sV(m*,)') + V ■ (7^V<) = inX 
«)'(^) = (/i)'(^) in dX 
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satisfy Adet[(V«:(x)),=i,2,3] = Ei=i,2,3 V«)'(a:) ■ [Vul+\x) x Vul+'{x)) > 

(we consider the expressions i + 1 and i + 2 modulo 3). We construct /] as the 
solutions of a system of ODE 

f = FH(/l).=l,2,3,^) 

fi\ £i 

J s\s=0 — Jo 

for an appropriate (-F*)j=i,2,3 defined below. 

We construct the functional F' : (C'=+2'°(5X))3 x [0, 1] C"^+2'°(9X) as 
follows. Assume X is a bounded domain. Let 70,7 G C'^+"^'^(X) and let 

7s = [(1 — s)7o + 57], s G [0, 1]. Assume < c < 70, 7 < C < 00. For s G [0, 1] and 
for (/Oi=i,2,3 G (C"=+2'°(9X))3 let be the C^-^^''^{X) solutions of 

V ■ (7s Vm*) =0 in X 
M*(x) = f^{x) in dX. 

Let be the C'=+2'"(X) solutions of 

'V- (7sVt;0 + V- ((7-7o)Vm^) = in X 
f*(x) = 5'*(x) in (9X, 

where G C'^'''^'"(9X) will be prescribed below. 

For this system let us consider the following auxiliary problems. Let A*, be the 
solutions of 

V ■ (7sVA^) = (Vu'+\x) X Vu'+\x)^ ■ V5i. in X 
A = in ax. 

From Theorem [321 A^^ G Lp{X) and {-fsdX'/du) G C^+2.°(aX) depend continu- 
ously on s. By integration by parts and summation we obtain 

J2 (Vm'+'(x) X Vm'+2(£) j . Vt;^(x) = J2 / ^'V ■ ((7 - 7o)Vm' 

1=1,2,3 i=l,2,3'^"^ 



Let (yf* = ^'jsdX'' /du, with /i chosen as 



E / A^V-((7-7o)V«0/ E / ( 

j=l,2,3'^"^ j=l,2,3'^^"^ 



Is r% 



■'+Hx 



Hence Ei=i,2,3 Vt^'(^) " (yu'+^{x) x Vu^ 

By defining -F*((/*)i=i.2,3, s) := (7* we can prove the boundedness and continuity 
of F' : {C^+'^''^{dX)f X [0, 1] ^ C"=+2'"(5X) as before. This yields a Theorem 
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analogous to Theorems 13.71 and 16.11 An important aspect for the argument to 
work is that we start with det[(Vu*(a;))j=i^2,3] > 1 and the evolution with 

maintains that property, hence J2i=i23ldx (^■5^) remains uniformly bounded 
away from zero. 

Theorem 6.2. Assume X is a (7'^+^'" bounded domain. Let 70,7 G C^+'^+^(X) 
and let 7^ = [(1 — 3)70 + S7],s G [0,1]. Assume < c < 70, 7 < C < 00. Let 
(/o)j=i,2,3 and 'Jq be chosen appropriately (e.g. 70 = 1 and /o(xi, 0:2, ^3) = Xi). 
Define F' : C''+^^''{dX) x [0,1] ^ C'=+2'°(5X) as above. Then, there exists a 
unique solution s ^ (/;) m C\[0, 1]; {C''+^''^{dX)Y) of the system of ODE 

§;fl = F'{{fs)i=l,2,3,s) 

£i\ £i 

Js\s=0 — Jo- 

For all s G [0, 1], the solutions u\ of (Pg) with corresponding boundary conditions 
fl are such that det[(VM*(x))j=i^2,3] > 1- 

This Theorem produces / = fs\s=i as the solution of the problem described at 
the beginning of this Subsection, with a condition involving finitely many equa- 
tions. 

Remark 6.3. In the definition of above, we could redefine /i = min(0,/i), 
resembling more closely the construction presented in Definition 13.31 Such a re- 
definition of provides a boundary condition {/*} with {{fl)'} of minimal L'^{dXY 
norm for each s G [0, 1], among all {/*} that produce non-decreasing determinants. 

Remark 6.4. The construction presented in this Subsection works in more gen- 
eral settings. We may consider X C M" and replace det[(V'U*(x))j=i^2,3] by 
I/[(Vti*(i;))j=i,...,m] for any multi-linear function L : (M")™ — > M. 

And if H : (M")'" — )■ M is a continuously differentiable function with differential 
DH{z) uniformly bounded away from zero in the set {z G (M")'" : H{z) > 1}, 
then the construction presented in this Subsection also works when we replace 
det[(VM^(x)),=i,2,3] by H{Vu%x),=i_m). 

Remark 6.5. Extensions for conditions involving multiple equations at finitely 
many points can also be addressed with this scheme. 
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